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Abstract. In this article, wc report on sixtii-ordcr and scvcntli-ordcr iterative mctiiods for solving nonlinear 
equations. In particular sixth-order derivative-based and derivative-free iterative families are constructed in 
such a way that they comprise a wide class of sixth-order methods which were developed in the past years. 
Weighting functions are introduced to enhance the algorithmic efficiency whereas an appropriate parametric 
combination gives weight-age flexibility in between those weighting functions. The usage of weighting factors 
and weighting functions define a wide class of iterative schemes for solving nonlinear equations. The freedom 
to construct different parametric combinations as well as different forms of weighting functions makes the 
iterative schemes more accurate and flexible, It means that one can easily modify the scheme by changing 
weight functions and parametric combination. 



1. Introduction 

Nonlinear algebraic equations are at the heart of many problems of nonlinear science. However, it is not always 
possible to find exact analytical solutions of these nonlinear equations which very often have to be solved numerically. 
Iterative methods are numerical algorithms which, starting from an initial guess in the neighborhood of a root of 
the nonlinear equation, manage to refine it and achieve convergence to the true root after several iterations . We 
can divide them in two classes, namely memory-based and without-memory iterative methods. From the point of 
view of efficiency and stability, multipoint iterative methods are superior than single-point methods. Kung and 
Traub conjectured [T3l that multipoint iterative methods without-memory have at most 2" order of convergence for 
n + 1 function evaluations per iteration. The further classification falls into derivative-free and non-derivative-free 
iterative methods. The multipoint iterative methods for nonlinear equations, which use derivatives, are efficient. 
The most famous derivative-based iterative method is Newton's method (NM) which is written as 

(1) ^n+l =x„ - , n = 0, 1, 2,3, ■ ■ • . 

It is quadratically convergent in some neighborhood of the root a of f{x) = 0. The computational effort invested 
in the derivative of a function is not always comparable to the evaluation of the function itself l2l l3l l4l [5l [6l 171 
|8]. Steffensen ^ gave a quadratically convergent derivative- free iterative method for nonlinear equations. The 
Stcffcnsen's method (SM) is given by 

!W„ = X„ - Kf{x„), 

Xn + 1 - X„ df (x„) ■ 

Clearly, df{x„) is a finite difference approximation of f'{xn). The SM approximation merely replaces the derivative 
evaluation by a function evaluation. Actually both iterative method NM and SM achieve the optimal order of 
convergence 2, which is in agreement with the conjecture of Kung and Traub. The computational efficiency of an 
iterative method (IM) is frequently calculated by the Ostrowski-Traub's 1141 1151 efficiency index 

(3) E{IM) = pi/'*, 

where p is order of convergence and d is function evaluations per iteration. The optimal computational efficiency is 
defined by expression 

(4) £;W=2"/("+l). 

We now briefly review a number of three step sixth-order convergent methods. Many three steps sixth-order 
convergent methods use three function and one derivative evaluations. Sharma and Guha [9] proposed the following 
sixth-order convergent family by introducing one parameter (SG) 



y?i — Xn 

(5) < Zn = Dn ■ 



/'(^,.)' 

/(^r.) /far.) 
f(xr,)-2S(y„) !'(x„)' 

f(x„) + af(yr,) f{z„) 
/{a:„) + {a-2)/(i/„) f'{x^)' 
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where a G SR. Neta [101 uni-parametric sixth-order family consists of three step (NTl) 



(6) 



/(a;„)+a/(t,„) J(iJ„) 



^ J/- 



/(x„) + (<i-2)/(!;„) 
/(=:«)-3/(!;„) /'(x„)' 



where a g SR. In 2011, Neta '11' provided a more efficient sixth-order scheme (NT2) 



(7) 



Vn 



Vn 
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'/(x„) /(x„) 



Chun and Ham |12l also presented a family of sixth-order iterative methods (CH) 

y-n. 



(8) 



^" /(^„)-2/(h„) /'(a;„)' 



where real valued function H{s) satisfies H{0) = 1, H'(0) = 2, and ^„ = f{yn)/f{xn)- Grau et al. 1161 developed 
the sixth-order (GR) 



(9) 



Vn 

Zn 



. „ _ fM f(y„.) 

y" /(3:„)-2/(h„) /'(a;„)' 



Alicia cordero et al. |17| constructed the following sixth-order scheme (AL) 



(10) 



Vn 


= X 


Zn 


= X 


Xn + 1 


= z, 



/lzn,!/nJ-|-/lzn,a;„,a;„J{z„-y„) ' 



where 9 g SR. Finally, Sanjay K. Khattri et al. |21| presented a paper. Unification of sixth-order iterative methods, 
which covers all the features of sixth-order iterative methods (SKI) 



(11) 



I 



aj I , where ai = 2, 



" - Pfej (l + ( "''""/IXo''"' ) ' -here b, = 2/^1 



where, {aj, fe^, /ii, ^2} 6 SR, ™ G Z+, and / G Z+ are independent parameters. Sanjay K. Khattri et al. [ISJ also 
developed four parameter sixth-order derivative free family (SK2) 

(12) 

Vn = X 



fM-f(x„-Kf(x„)) 



1 + 



/(I/,.) 



/{a:„-K/(a;„)) 

f(x„-Kf(x„)) 



f(xn-nf{x„)) 

f(y-n) 

f(x„-iif(x„)) 



' S(y^) 



In 1191 . R. Thukral constructed Steffensen-Six.one method (TSl) and Stcffensen-six.two method (TS2), which are 
respectively 



(13) 



Wn =Xn+f(Xn) 
Vn — Xfi 



/{™„)-/(x„) h 



Xn + l 



1 - A 



-1 nihil 



fM-Hy^) 
f{Z]-fly^) 



fiVn), 
/(2n), 
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and 



(14) 



x„+i 



■ Xn + f(Xn), 



■ Vn 



1 - A 



/(Zr.)-/(H„) 



/(S/n), 



where A G 5fi/{0}. F. Soleymani et al. 20 discussed the following sixth-order derivative-free method (FSl) 



(15) 



Vn 

Xn+1 



■ X„ + f{x„), 



1 1 f(Vn) I 1 I o f(Vn) 
, l-|-/[a:„,m,.] /(z„) 



where f[yn,Zn] = •^^^^j-yl^"'' ' /[^'n'"""] 
derivative-free scheme (FS2) 



J(ij),.)-/(3:„) 



F. Soleymani [22] also presented an other sixth-order 



(16) 





— Xn 


Vn 


— Xn 




= yn - 




— 



/[™„,2„]-|-/[z„,!/„]-/[m„,y„] ' 



where f} G 5(i/{0}.In 2013, F. Soleymani I26| developed some efficient seventh-order derivative-free methods. First 
is (FS3-1) 



(17) 





= Xn-\- 


Vn 


— Xn 




= yn- 


Xn 


— Zji 



flXn,yn]+flVn,Wn]-f[Xn,W„] ' 



where 7 and 5 are real numbers. Second is (FS3-2) 



(18) 





— Xn 




— Xn 




= yn- 


Xn 


— Zji 



f{Xn), 



f(y^) 



f[Xn,y-n] + f[yTi,,Wn]- f\xT^,Wn] ' 



-I , f{y^) , fU^) 



where p and r are real numbers. Third is (FS4-1) 
(19) 

On =Xn + f{Xn), 
f(.x„) 



Wn 


= Xn + 


Vn 


— Xn 




= yn - 


Xn 


— Zn 



/(Vn.) 

flXnyynl + flynylUnl—flXntWn 



^ f{^n) ^ f(y^) ^ 1 (i-l-/b„,»„])^ f I f{x^) ] '^if(x^) ^° f(vj^) 



2-f[x„,-w„] 
(l-flx„,w„]y^ 



f(x„) 



'f(x„) 



/[uj„,z„l 

where lo and </> are real numbers. Fourth is (FS4-2) 



+ 7{f^ + + {2 + /I-"' -"](3 + flXn,Wn])} j 



' f(Xn) 



(20) 



Wn 


— Xn 


Vn 


— Xn 




= Vn 


Xn 


— 



f(x^) 



" ~ ^" flx„,y„]+f[y„,m„]-f[x„,w„] 



^ + yfer + jfery + + f[xn,wn]{-3 + /k, «;„])} 



f(yn) 
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2. A More general formulation of sixth-order iterative methods 

In this section, we go beyond the iterative method (SKI) and present a more general formulation which covers 
most of the sixth-order iterative methods based on one derivative evaluation f'{x). The proposed iterative scheme 
is (FDl) 



(21) =yn-G{ti)j}^, 



where ti = , t2 = j^^"? and ts = {^^"j . we state the following theorem. 

Theorem 1. Lot / iDCSR— ^JRbea sufficiently differentiable function, and a € D is a simple root of fix) = 0, 
for an open interval D. If xq is chosen sufficiently close to a, then iterative scheme II21II converges to a. If G and 
H satisfy 

(22) G(0) = 1, dO/dtiijd) = 2, /i'{0,0,0) = 1, dH/dti(Q,Q,Q) = 2, 

and both G, H have bounded higher order derivatives then the iterative scheme II21II has convergence order at-least 
six and along with conditions II22II if 

(23) 9^/^/9*? (0,0,0) = ePG/dtlijd) + 2, 9///9t3 (0, 0, 0) = 1 
then the iterative scheme 12111 has convergence order at least seven. 



Proof, let a be the simple root of f{x) = and e„ = Xn — a. Further we denote ci = f'{a) and c^. = f'(a) ' 
where fc = 2, 3, ■ • ■ . The Tyalor's expression of / and /' around a in terms of e„: 

(24) f{xn) = ci [en + C2el + cge^ + c^el + c^e^ + c^el + 0(el)] 

(25) f'{xn) = ci[l -h 2c2e„ + Zc^el + Ac^el + Scge^ + ecge^ + 0(efj] 

Substituting 1241 1 and 1251 1 into y„, we obtain 

_ fi^n) _ f{Xn) _ _ f{Xn) 



/'(^n)' f'(x„) f'iXn) 

■^)el + (3C4 - 7C2C3 + 4ci)e* + (4c5 - 1 
(26) -I- (-17C3C4 + 33c2C§ - 52c3ci -I- 28c4C^ - 13c2C5 -|- 5c6 -|- 16c|)e^ + O(eJ^) 



C2el + (2c3 - 24)el + (3c4 - 7c2C3 -I- 4ci)e* + (4c5 - IOC2C4 - 6c§ + 20c3ci - 8ct)e^ 



Expansion of f{yn) around a gives 

fiVn) = ci[c2el - 2(-C3 + cl)e\ + (3c4 - 7c2C3 + ^cl)e\ - 2{-2cs + 5C2C4 + 3c| - 12c3c2 + %c\)el 
(27) ci(37c2ci - 73c3ci + 28c| + 34c4ci - I7C3C4 - I3C2C5 -I- 5c6)e^ -|- 0{el)] 

We define G(ti) = 1 + 2ti + Mot^ + Mitf + M2tf + higher powers in ti , where Mq, Mi and M2 are finite real 
numbers. By using 124II - 127II . the expressions for Zn and f{z 



Zn — Oi = —C 



= -C2(-5ci + Mo4 + C3)ei + (-2C2C4 - 2ci + 32c3ci - 36ct - Mi4 - 6Mo4c:i + lOMo4)el 
(-12MoC2c| - 9Mocic4 -I- 74Moc|c3 - 62Mo4 - 7c3C4 -I- 66c2C§ - 262c3c| -I- 48c4c| - 3c2C5 

(28) + 170c| - M24 + 13Mic| - 8Micic3)e^ + O(e^) 

f{zn) = -ciC2(-5ci + Mo4 + C3)ei + ci(-2c2C4 - 2ci -f 32c3ci - 36c^ - Mi4 ~ 6Moc|c3 

+ 10Moc^)e^ -I- ci(-12M'oC2C§ - 9Mocic4 -I- 74Moc|c3 - 62Mo4 - 7c3C4 + 66c2C§ - 262c3c| 

(29) + 48c4ci - 3c2C5 + 1704 - M24 + 13A/ic| - 8Micic3)efj -I- O(eJ^) 
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Now, we define H{t-i,t2,t-^) = 1 + 2ti + Rot2 + R\t-j, + R2t\ + Rzt\ + R^t^ + higher powers in ti, t2, t^. Again, 
substituting |[25ll . |[28t and |[29t in H2TI . we obtain 

x„+i - a = Zn-a- H(ti,t2,t3)— 

= -C2{-hc% + Moci + C3)(HiMoci - 5Rici + 6ci - R2c\ + Rxcz - C3)e® 

+ (IOR0C2M0 - Rf)C%Ml + lOiJiMicI - 172i?ic^Mo + 20iJiM^c| + 2Mocic4 + 6Mocic§ 

- 68M0C2C3 + 88Moc^ + 124:RiMo4c3 - URiMoc^cl - WRiM^c^cs - 2RiMi4Mo - 2RiMi4c3 

- 2i?oc|MoC3 - iRiaMocl - 22c4c| - 66c|c^ + 366c3cf + 2c| - 356c| + 4040203 - 6Mic| 

- 25i?6C2 + lOiJeCaMo + 10R6ctc3 - flec^M^ - i?6C^c§ + Mic|c3 + 76ii;2C^ - 2R6c|MoC3 

- 4c2i?iC4C3 - 342R1C34 + 20ii:iC4c| + 10fioctc3 + 64c|i?ic§ - Roc|c§ - 25i?oc^ + 360ii;ic^ 

(30) - 2iJic| - 60i?2c|c3 + 6i?2c|c§ + 2R2cic4 - 18iJ2c|Mo + i?2C^Mi + 10i?2ctMoC3)e^ + 0(6^) 
By equating one of the factor of coefficient of to zero, we have 

(31) RiMoci - 5Rici + 6c^ - R20^ + iJiC3 - C3 = 0. 



Now by equating the coefficient of 03 and Cj in II31I I. we find 

Rl — 1 = 0, (coefficient of C3) 



(32) 



RlMo - 5Ri + 6 - i?2 = 0, (coefiicient of c|) 



where Mq = 1/2 cPG/dtl{0), Ri = dH/dt3{0, 0, 0) and R2 = 1/2 d'^H/dtl{0, 0, 0). Which completes the proof. □ 

Remark 2.1. It is obvious that sixth-order family II21I I can construct all the presented iterative methods in l(5 ll -l|9 ll 
and ((TTJ. 



JTO} can be written as (FD2) 
(33) 



Vn 

— X 



where ti = and A(ii) = (^^j^^ , with 9 G R and /[z„, = ^'^"'^"i/'^'^"'' ■ The more general form 

of II33II is stated in the following theorem. 

Theorem 2. Let / iDCSR— >3(ibea sufficiently differentiable function, and a G D is a simple root of f{x) = 0, 
for an open interval D. If xq is chosen sufficiently close to a, then iterative scheme II33I I converges to a. If A satisfies 

(34) ^(0) = 1 

then the iterative scheme I I33I I has convergence order at least six and if along with I I34I I. it satisfies 

(35) dA/dti{(}) = 2 
then 13311 has convergence order at least seven. 

Proof. We define A{ti) = 1 + Moh + Mitf + M2tf + higher powers in s. By using the 12411 - 11271 1. we obtain 
Zn-a = yn-a-{l + Moti + Mi4 + Mzt?)^^^ 

f (Xn) 

= {-ol(-2 + Mo))el + C2(-9c| -|- 7Moc| - Mio?, + 703 - 4MoC3)e^ -|- (IOC2C4 -|- 6c| - 44c3ci 
-I- 30c| - 6M0C2C4 -I- 38Moc|c3 - 33Moc| - M2C^ - 6Mic|c3 -I- 10Mic| - 4Moc§)e5 -|- (17C3C4 

- 70c2c| + 188c3C^ - 62c4c| -I- I3C2C5 - 88c| -f 74Mic^C3 - 9A/ic|c4 - 12AfiC2c| - 62Mic^ 

+ 13Af2ci - 8M2cic3 + 129Moc| - 12AfoC3C4 + 68A/oC2C§ - 225AfoC3c| -I- 55AfoC4C^ - 8AfoC2C5)e® 

(36) +0{el), 

f{zn) = ci[(-ci(-2 + Mo))el + C2(-9ci + 7Mocl - Mic^ + 7c3 - 4MoC3)ei + (IOC2C4 + dcj - 44c3ci 
-I- 30c^ - 6M0C2C4 -I- 38Afocic3 - 33Afoc2 - Af2ct - 6Aficic3 -I- 10Afic| - 4Moc|)e^ -I- (I7C3C4 

- 70c2C§ + 188c3ci - 62c4c| + 13c2C5 - 88c| + 74Aficic3 - 9Afic|c4 - 12AfiC2C§ - 62A/ic| 

+ 13Af2ci - 8Af2c|c3 -I- 129Afoci - 12AfoC3C4 -I- 68AfoC2C§ - 225Afoc3c| -I- 55AfoC4c| - 8Afoc2C5)e® 

(37) +Oiel)]. 
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By substituting values from I I24I1 - II27I I. I I36I I and I I37I I in x„-^-l from II33I I. we get 

fiZn) 

- « = z„ - a - — — — - 

e„+i = ci(-2 + Mo){Moc^ + 2c3 - 2ci)e^ - ci(364 - 46Moct + 4Mic^ - 2MictMo + liM^4. 
(38) - 10M,?C3c| - 2Mic|c3 - 64c3c| + 57Moc|c3 + 6C2C4 - 3MoC2C4 + 22c§ - 12Moc§)eJ^ + O(e^). 

Clearly, if dA{t\) / dt\{Q) = Mo = 2 then convergence order is seven otherwise is six. □ 




Now we present derivative-free optimal fourth-order iterative schemes (FD3): 

— Xti K, f iXji') ^ 

(39) {V" =Xn- > 

= y-- [7I^^0(*l) + 7I^^l(*2) + 7pf^G2(ti,t2)] , 
where t\ = t2 = J^^^- The convergence of II39I I is given in the following theorem. 

Theorem 3. Let / :Z)C5R— >3tbea sufficiently differentiable function, and a G -D is a simple root of f{x) = 0, 
for an open interval D. If xq is chosen sufficiently close to a, then iterative scheme II39II converges to a. If Go(ti), 
Gi(t2), G2{ti,t2), go, gi and 92 satisfy 

r Go(0)) = f, Gi(0) = 1, G2(0,0) = 1, 

(40) I dGo/dti{0) = 1, dGi/dt2iO) = 1, 9G2/9ii(0,0) = 1, 90/9*2(0,0) = 1, 
[ go + gi + 92 = 1 

and all higher order derivatives of Gi, i = 0, 1, 2 are bounded then the iterative scheme I I39I I has convergence order 
at least four. 

Proof. Error term at nth-step is denoted by e„ = Xn — a. we call ci = f'{a) and Ck = ^ f'(j^^ ' where A; = 2, 3, ■ ■ ■ . 
Taylor's expressions give: 

(41) f{xn) = ci[e„ + C2el + c^el + C4e^ + cse^ -|- cee^ +0(e^)]. 

/(w„) = — ci(-l -I- KCi)e„ -I- ciC2(— 3kci -|- 1 + K^cf)e^ - ci(4kciC3 — C3 - 3c3K^Ci -|- C3K^cf -|- 2c2KCi 

— 2K'^cfc|)e^ + ci(— 5kciC4 — 5c2KCiC3 + 8k^c^C2C3 + K^cfc^ + C4 -|- 6c4K^cf — 4c4K^c'1 + c^k^c^ 

(42) - 3c2K3cfc3)e4 +■■■ + O(e^). 

By using Mil l and II42I I. we get 

f(x„) — f{Wn) 

f[Xn,Wn] = — = -ClC2(-2 + KCl)e„ -|- Cl (3C3 - 3kCiC3 + C-J,K?c{ — C2KCl)e^ — Cl( — 4C4 

Xn 

pn '^3 2 2 \3 / 2 2 2 2 2 

+ 4c2ftciC3 -I- 6kciC4 — 4c4K Cj + C4K cf — 2k CjC2C3)e„ -|- ci(2k Cj^C3 -|- 8k Cj^C4C2 

— IOKC1C5 -f 5c5 — 7c4KCiC2 — 3KC1C3 -I- IOC5K Cj — 5C5K cj -|- C5K Cl — 3c4K CJC2 -|- k CjC2C3)e„ 

(43) + ■ ■ ■ + 0{el) 

By substituting II41II and II43II in y = x — f{x„)/f[x„,'Wn], we obtain 

yn — a = — C2( — 1 + KCi)e5; + (2c3 — 3kciC3 + czk^c\ + 2c|kci — 2c| — K^c\c^)e^ + (3c4 + IOC2KC1C3 — 6KC1C4 

(44) -I- icAK^c\ - CAK^cl - 7k^CiC2C3 - 7c2C3 - SkciCj -I- 2c2K^cf C3 -|- 3K^cf -I- 4c2 - K^cf C2)e^ H h O(e^) 

Taylor's Expansion of / around ?/„ by using 1144 l l is 

fiVn) = ci(-C2(-l + KCi)e^ -I- (2c3 — 3kciC3 + C3K^cf -I- 2c|kci - 2c2 — K^cfc|)e^ -I- (3c4 -I- IOC2KC1C3 - 6KC1C4 

(45) + ic4K^c1 - C4K^cf - 7K^cf C2C3 - 7c2C3 - 5kcic| -|- 2c2K^cf C3 -|- Zk?c\c\ + ic^ - K^cf c|)e* H h 0(ej^)) 

Now, we define 

(Go(ti) = 1 + L\t\ + higher powers in ti, 

Gi(i2) = l-f t2 + -I- higher powers in t2, 

G2(ti,t2) =l + ti+t2+ Nit^ + N2tl + higher powers in ti, t2, 
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where Li, Mi, Ni and N2 are real numbers. By substituting II41I I-I I45I I in x„+i = Vn — ^ f [ y^°m ] Go(ti) + 
7[^Gi(t2) + j^-S^G2{h,t2)] , we obtain 

En+i = (-1 + SO + 91 + ff2)c2{-l + KCi)el + 

(-1 + 90 + 31 + 32)(-2c3 + 3KC1C3 - C3K^Ci - 2c|kci + 2cl + K^cjcl)e^ + 

(-3goC4 - SQci + 6goC2C3 - 7K^cfc2C3 - 6KC1C4 + 3c4 + 892^^1 ci^ci - 3g2A^iC2K^Ci 

+ 52A^iC2K^cf + g2N2C2KCi + giMic^KCi + goK^cfiicI + 3go-f^iC2Kci - 3goLiC2K^cf + 3K^cfc2 

+ 4C4K^cl - C4K^Ci - 5KC1C2 - K^cf C2 + IOC2KC1C3 + 2C2K^cf C3 - 7C2C3 + 4C2 - 5S2KC1C2 + giK^cf Cj 

+ 6goKCiC4 - 4goC4K^Ci + goc4K^cf - 2goKciC2 + 2goK?c\c2 - g2N2C2 - giMic^ - g2NiC2 - goLic^ 
— 8goC2KCiC3 + 6go«;'^cf C2C3 — 2goC2K''cfc3 + 6giKciC4 + 6giK^CiC2C3 — 8giC2KCiC3 — 2giC2K^cJc3 
+ 6g2K^cf C2C3 - 8g2C2KciC3 - 2g2C2K^cf C3 + 6giC2C3 + 6g2C2C3 - 3giC4 - gic^ - 3g2C4 + 52C2 
(47) — giK^cf C2 — 4giC4K^Ci + giC4K'^cl + 6g2KCiC4 + 3g2K^cf — 4g2C4K^cf + g2C4K^cf )e* + ■ • ■ + O(eJ^) 

Clearly, if go + Si + 92 = 1 then the coefficients of and ej^ are zero and hence the convergence order is four which 
is optimal in the sense of Kung and Traub. □ 

Remark 2.2. In II39I I. if we equate go = Oi 9i = 0, 32 = 1 and define G2{ti,t2) = 1 + ti + £2 + «tf + /9t2' then 
the resulted iterative scheme is the first two steps of iterative method II12I I. Similarly if we equate go = 0, g2 = 0, 
gi = 1, K = —1 and define Gi(t2) = (1 ~ -^~^t2)~'*'i then it reproduces the ffist three step of iterative schemes I ll3t 
and |[T4t . 

In order to extend the iterative scheme II39I I. we add one more step and get (FD4) 



(48) 



Zn = y- 



5o (i 3 , i4 , 1 5 ) + 7[i^ 5l (t 1 , 1 3 , 1 4 , 1 5 ) + S2 (t 2 , t3 , t4 , i 5 ) 

+ 7P^^5'3(i2, i3, t4, ts) + j[^^^S4(il, is, t4, ts) + /[^|^%„] -^S (tl , fa . ^3 . ^4 , ts)] /(^n), 



where ti = j^^, h = f(ZZ\ ' *3 = T^f^, *4 = J^^*^') , *5 = and hi, ^2, /i3. ^14, hs. « 7^ are real numbers. 

The convergence of II48II is discussed in the following theorem. 

Theorem 4. By including the statement of (Theorem 3), iterative scheme 1148 II has convergence order at least six 
if 



(49) 



5o(0,0, 0) = 1, Si(0, 0,0,0) = 1, 52(0,0,0,0) = !, 53(0,0,0,0) = !, 54(0,0,0,0) = !, 
55(0,0,0,0,0) = !, 95i/9ti(0,0, 0,0) = 1, 952/9*2(0,0,0,0) = !, 
953/9*2(0,0,0,0) = !, 954/9*1(0,0,0,0) = !, 
955 /9*i (0, 0, 0,0,0) = !, 955 /9*2 (0, 0, 0, 0,0) = ! 



and all higher order derivatives of all functions are bounded. 

Proof. Taylor's expansion of various orders for 5;, i = 0, !,••■, 5 is given below. 
(50) 

'5o(*3,*4, *5) = ! + ai*3 + 02*4 + a3*5 + higher powers in *3, *4, *5, 

5i(*i,*3, *4, *5) =! + *!+ bitf + 62*3 + ^3*4 + t'4*5 + higher powers in *i, *3, *4, *5, 
52(*l,*3, *4, *5) = ! + *2 + 65*2 + ''6*3 + ^7*4 + bgts + higher powers in *2, *3, *4, *5, 
53(*2,*3, *4, *5) = ! + *2 + feg*2 + 6l0*3 + bii*4 + bi2*5 + higher powers in *2, *3, *4, *5, 
54(*i,*3, *4, *5) =! + *!+ fei3*i + fei4*3 + ^15*4 + tii6*5 + higher powers in *i, *3, *4, *5, 
'S'5(*l,*2, *3, *4, *5) = ! + *i + *2 + fci7*f + ''IS*! + ''19*3 + ''20*4 + ''21*5 + higher powers in *i, *2, *3, *4, *5, 

where ai, a2, ag, bi, 62, 63, ''4, ''5, ''6, h, bg, bg, bio, hi, f'12, 613, ''14, ^is, ''16, ''17, ^is, ''19, ''20, and 621 
are real numbers. By using calculations from II4!| |-| |47| | for x^+i = z„ — (i-fei-ft2-fa3-fa4-^5) ^^(^f^^f^^ ^g) _|_ 

j[^^i^5l(*l,*3,*4,*5) 

+ j. ; ,^^,„] 5'2(*2,*3,*4,*5)+ ; ; ,^%„ ]g3(*2,*3,*4,*5)+ j.[^^"^^] g4(*l,*3,*4,*5)+ j[,^^^„] 55(^ 
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we obtain tlie following error equation. 

O 

+ ^5b2l92C2KCi — 2hr:,b2lLiC2KC\ + /l562lilC2«; Cj^ + 5c2/i5 + hiK C1C2 + /14KC1C3 — 4/l4C2«;Ci + hiK C1C2 



+ hlKClC3 — 3/iic|kCi — c|h2tCl — 2C2/13KC1 — bc^h^KCl + /l2KClC3 + h^KClC^ + h^K?c\c2 + /l5KClC3 

— a^K?c\c\ — a^KCic^ + 4a3C2KCi + 2a3hi(;2C2 + ■ ■ ■ + /!.3f)i2/^ic|/t^Cj — 2/i2bg(?ic|KCi + h2bg,g\K? c\c^ 
+ '*2''892C2KCi — 2/i2b8Lic|KCi + /i2b8-^^iC2's^Ci)c2(— 1 + KCi)(— giK^cJc^ — Lic|k^Cj — g2N-i<^K? c\ 



+ K''cfc2 + giLiC2K''cf + g2LiC2K'^cl — 4c2KCi + 2giC2KCi + 2g2-'ViC2KCi — 32C2KC1 — 2giLiC2KCi 
+ KC1C3 + 2Lic|kci - 2g2iiC2KCi - C3 - g2A^2ci - g2iVici + 2g2ci + 3c| - giMiCj - Lic| + giLic^ 
(51) +S2i^ici)e6 +0(e^) 



□ 



Remark 2.3. In II48I I. by equating ft = —1, gi = 1, 92 = 0, = 1, hi = h2 = = = and defining 
G2(*2) = (1 — A~lt2)~^, 53(*2, *3, t4, is) = (1 — A~ 1*2) gives the iterative scheme llT3t . by setting k = — 1, gi = 1, 
g2 = 0, hi = h2 = = h4 = = and defining G2(t2) = (l~A~-'^i2)~'^ ,5o(t3, t4, is) = 1, we obtain the iterative 
schemes II14I I. II12I I can be constructed by choosing parameters gi = 0, 92 = 1, /is = 1, /ii = /i2 = /13 = '14 = and 
weight functions G2(ii,i2) = 1 + ti + t2 + +134, S'5(ii , i2, is, t4, is) = 1 + ti + i2 + +13^ +^*5- 

F. Soleymani constructs derivative-free sixth-order iterative scheme III6I I. In the following theorem, we present a 
modified iterative family with convergence order seven. 

Theorem 5. Lot / iDCJf— >Jffbea sufficiently differentiable function, and a g D is a simple root of f{x) = 0, 
for an open interval D. If xq is chosen sufficiently close to a, then iterative scheme (FD5) 



(52) 



where 







Vn 


— 




= Vn 




— 



H(t3,t4, is) f[^^^y^^_(^h-l)flz„,m„] + hflJ'„"J„]-hfly„,x„] + {h-l)f[y„,m„] 



fiVn), 



(53) 



1 + ii + aitj + higher powers in ti, 



1 + t2 + a2t2+ higher powers in t2. 



Go(ti) 

Gl(i2) 

G2(ii,i2) = 1 -f ii + t2 -f a3ii + i4i| + asili2 + higher powers in ii, t2, 
H{t-j, t4, is) = 1 + a6i3 + ayt-i + a8is + higher powers in tg, t^, is, 

1, 2, .., 5 are consistent with previous definition in II48I I. gi, 92, aj , j = 1, ■ ■ ■ ,fi are real numbers. 



(54) 



has convergence order at-least six if 

JGo(0) = 1, Gi(0) = l, G2(0,0) = l, //(O, 0, 0) = 1, 
\(iGoMi(0) = 1, dGi/dt2{0) = 1, 9G2/9ii(0,0) = 1, 9G2/9i2(0,0) 

and has convergence order seven if along with II 54 I I. the following condition is valid 

(55) 9/^/9*5(0,0,0) = 0. 



Proof. The terminology of (Theorem 3) and expressions in I I4H . II42I I assist to perform all necessary calculations by 
using Maple 13 I23| software-package in order to get the following error equation: 

e„_|_i = a8C2(— 1 + k,ci)(k? c\c2 — C3 -f KC1C3 -|- 3c2 — 4c|kci -|- giaic^ — aiC2+ 2aic'^KCi — aic|K^cJ — (/2C2tci 

+ 2giC2KCi - giK^CjCj -I- 32aic| - gia2C2 - g2a3C2 - 32asci - 92a4ci -I- 2g2C2 - 2giaiC2KCi + giaic^K^c^ 

- 2g2aic^KCi + g2aiC2K^cl + 2g2a-iC2KCi ~ g2a'iC%K^c\ + g2azc%K.ci)'^ 

/ 

H (k C1C2 — C3 -I- KC1C3 -1- 3c2 — 4c2«;ci + giaiC2 — aiC2 -I- 2aiC2KCi — aiC2K Cj — g2C2K.ci + 2giC2KCi 

- giK^cf c| + g2aiC2 - g\a2C2 - g2as,c% - g2asC2 ~ 92a4c| + 2g2C% - 2giaic%KCi + giaic%K'^c\ 



A A. A ^ O O O A 

2a8C3 — 30a8C2 + 3aeC2 — Za^K C1C2C3 + agcsK C1C2 + 3aBKCiC2C3 — baegiK c^C2 
4 I A „3„3A I A 2„2„4 



(56) 



- aeaiC2 
+ 2ae,giKCiC2 
+ 4:a6ai4K^cl)el^ + 0(efj 



AaagiK'' C1C2 + Aafig2K,^ C1C2 — ba(,g2KCiC2 — a6(/2K^cfc2 — aggiK^^c^ — aaai^K.'^cf 



IMPROVING THREE-POINT ITERATIVE METHODS FOR. SOLVING NONLINEAR EQUATIONS 
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Clearly, ag = dH/dt^(0,Q,0) is a factor in asymptotic error constant of ej^, if ag = then I I53I I has convergence 
order seven otherwise convergence order is six. □ 

Remark 2.4. Iterative scheme II53I I reproduces F. Soleymani derivative-free sixth-order iterative scheme II16I I if 
(57) = 0, 91 = 92 = 0, Go(ii) = 1, H{t3,U, tg) = 1. 

The derivative- free version of II33I I is (FD6): 



(58) 



where 



(59) 







Vn 


— ^'n 




= Vn 




— Zfl 



/[2n,H„] + (/[z„,a:„l-/[ii.„ 



Go(ti) 

Gl(t2) 
G2(tl,t2) 



1 + ti + ajtf + higher powers in ti, 
+2 . 



1 + tl + t2 + a^tl + a^tl 



+ ^5*1*2 + higher powers in ti, t2, 
1 + agts + a^t/^ + agts -I- higher powers in ta, i4, is, 



ii, i = 1, 2, ■ • • ,5 are consistent with previous definition in II48I I. gi, (/2, h, aj , j = I, ■ ■ ■ ,8 are real numbers. 

Theorem 6. Let / :DC5R— >5Rbea sufficiently differentiable function, and a g D is a simple root of f{x) = 0, 
for an open interval D. If xq is chosen sufficiently close to a, then iterative scheme I I59II has convergence order 
at-least six if 



(60) 



100(0) = 1, Gi(0) = l, G2(0,0) = l, /^(O, 0, 0) = 1, 

\dGo/dti(0) = 1, dGi/dt2(0) = 1, 9G2/9ii(0,0) = 1, 902/9*2(0,0) 



1. 



Proof. We obtain the following error equation 
e„+i = (-1 -I- kci)[k?c\c% - c-j, + KCiC'i + 3c| - 4c2KCi + giaiCj -I- 920.10% - gia2c| - g2C%a^ - g2a4C2 

- 920302 - aicj - 2giaic%KCi + giaic^K^cl - 2g2aic%K.ci + 9201 clre^cf -|- 2g2a3ciKCi -|- 2aic|KCi 

— a\C2K, Cj -f 2(/iC2KCi — giK C]^C2 — g2C2KCi + 2g2C2 — g2a3C2K c-^ + g2C2a^K.ci){K c^C2 — C2KC1 + 3c2a8 

- ages + 2c|agg2 - aic|ag - c|aga2gi — c|a8g2a3 — c%ag,g2a5 — c|agg2a4 -I- ciagcsK + aic|aggi 

+ aiC2agg2 — cfc|agK^gi -f- 2c\c\asK,gi — c\c^Kag,g2 — c\c\k^ ag,g2a3 + ciC2Kagg2a5 + 2ciC2«;agg2a3 

— aiC]^C2agK + 2aiciC2agK + aiCiC2a8K gi + aiCiC2K asg2 — 2aiciC2asKgi — 2aiciC2Kasg2 + CiC2asK 



(61) 



■ 4ciC2ag/t)c2e^ + O(eJ^) 



□ 



3. Numerical computations 



Definition 1. Let a;^— 1, Xn and Xn+i be successive iterations closer to the root a of f{x) = 0, the computational 
order of convergence (COC) | 24| . can be approximated by 

ln\(xn+i - a){xn - «)~-'-| 



(62) 



COC ! 



ln\{x„ - a){x„-i - a) l| 



A set of twelve functions is listed in Table 1 which is taken from I25| to validate the iterative methods and 
performance. In all methods, twelve total number of function evaluations (TNFE) are used. 



3.1. Sixth-order convergence performance evaluation of FDl. For the purpose of comparison between FDl 
and SG, NTl, NT2, CH, GR, AL, wc derive an iterative method from (I2TJ which is given as (FDl-Ml ) 



(63) 



where ti = {fr"? , t2 




Zr, - Hi{ti)H2{t2)H3{t3)j^, 

= 7-57-177. ■f^2(*2) 



1 



1 -I- 2*2 and H3(t3) 



1 



l-2ti '--i^--^' l_2ti-t^ 

1—1 itg • Table 2 shows the performance of different listed iterative methods in terms of absolute error (ja;„ — a|). 
Clearly FDl-Ml is competitive with other developed methods and freedom to choose different weight functions, 
makes it more accurate and stable in the third step of iterative scheme II21I I. 
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Table 1. List of test functions 



Functions 


Roots 


fi{x) = exp{x) sin{x) + ln{l + x'') 


« = 


f2{x) = + + 4x2 _ 15 


a = 1.148538... 


/3(a;) = (x - 2){x^° + x + 1) exp(-x - 1) 


a = 2 


Si{x) = exp{—x^ + X + 2) — cos{x + 1) + a;^ + 1 


a = -1 


fs^x) = (x + 1) exp{sin{x)) — x"^ exp{cos{x)) — 1 


a = 


f%{x) = sin(x)'^ — x'^ + 1 


a = 1.40449165... 


frix) = 10 exp{—x^) — 1 


a = 1.517427... 


fs{x) = (x2 - - 1 


a = 1.414214... 


fg{x) = ln{x'^ + X + 2) -x + 1 


a = 4.15259074... 


fio{x) = cos(z)2 - x/5 


a = 1.08598268... 


fli{x) = a;!" - 2x3 - X + 1 


a = 0.591448093... 


fl2{x) = exp{sin{x)) — x + 1 


a = 2.63066415... 



Table 2. Numerical comparison between absolute errors(|x„ — a\,TN FE = 12) 



/„(x),xo 


FDl-Ml 


SG 


NTl 


NT2 


CH 


GR 


AL 






a = -1 


a = -1 








e = -1.01 


fi, 0.25 


2.71 e-142 


2.82e-121 


1.22e-127 


4.25e-84 


5.12e-82 


9.23 e-93 


6.22e-141 


/2 , 1.1 


3.90 e-190 


3.72 e-153 


1.23 e-122 


1.47 e-95 


2.45 e-104 


4.43 e-141 


1.60 e-123 


/3, 2.1 


2.31e-181 


3.75e-135 


6.27e-135 


2.89e-93 


1.07e-94 


6.52e-109 


1.80e-147 


/4, -0.5 


5.28 e-218 


8.50 e-181 


5.30 e-143 


2.70 e-173 


1.27 e-174 


5.89 e-207 


8.99 e-191 


/s, 0.25 


1.96e-247 


4.05e-184 


2.51e-165 


3.77e-213 


7.33e-162 


2.97e-243 


2.20e-215 


/6, 1-2 


3.79e-215 


1.41e-181 


7.22e-173 


3.90e-123 


1.69e-123 


6.80e-154 


2.42e-180 


/7, 1.0 


5.26e-123 


2.88e-95 


2.10e-113 


3.87e-69 


6.28e-67 


1.81e-76 


2.55e-118 


fs, 1-6 


2.27e-181 


1.18e-157 


1.65e-97 


6.69e-59 


1.23e-103 


2.26e-147 


4.14e-60 


/9, 4.4 


2.92e-402 


3.136-366 


4.80e-387 


2.35e-344 


6.52e-341 


1.89e-353 


2.41e-432 


flO, 1-5 


2.33e-56 


2.22e-53 


3.42e-64 


3.58e-31 


10.1 


5.45e-56 


8.42e-58 


fii, 0.25 


8.99e-115 


1.30e-97 


1.54e-91 


1.31e-53 


3.86e-45 


4.77e-89 


1.72e-92 


/i2, 2.0 


3.53e-173 


1.80e-132 


3.42e-138 


1.84e-117 


2.78e-115 


3.52e-126 


8.05e-162 



Table 3. Computational order of convergence (COG) 



fn(x) 


FDl-Ml 


SG 


NTl 


NT2 


CH 


GR 


AL 


h 


5.9999 


5.9994 


6.0003 


5.9989 


5.9987 


5.9995 


5.9965 


f2 


6.0000 


6.0000 


5.9999 


5.9997 


5.9998 


6.0000 


5.9933 


fs 


6.0000 


6.0000 


5.9999 


5.9995 


5.9996 


5.9999 


6.0012 


u 


6.0000 


6.0000 


6.0000 


6.0000 


6.0000 


6.0000 


5.9998 


k 


6.0000 


6.0000 


5.9999 


6.0000 


5.9999 


6.0000 


5.9999 


h 


6.0000 


6.0000 


5.9999 


5.9999 


5.9999 


6.0000 


5.9998 


h 


5.9999 


5.9992 


6.0001 


5.9974 


5.9972 


5.9991 


5.9860 


fs 


6.0000 


6.0000 


6.0001 


5.9985 


5.9999 


6.0000 


6.0971 


/9 


6.0000 


6.0000 


6.0000 


6.0000 


6.0000 


6.0000 


6.0000 


/lO 


5.9806 


6.0038 


6.0098 


5.8343 


3.4125 


5.9880 


6.1207 


fii 


5.9999 


5.9981 


6.0001 


5.9909 


5.9856 


6.0000 


6.0120 


fl2 


6.0000 


6.0000 


6.0000 


6.0002 


6.0002 


6.0001 


5.9981 



3.2. Seventh-order convergence performance evaluation of FDl. To compare the results, we construct 
seventh-order iterative scheme (FD1-M2) from II21I I: 



Vn — Xfi 



(64) {z„ =yn- G(ti) 



r„+i =z„-Hi{ti)H2{t2)H3{t3)j 



where t'-s are defined previously, G(ti) = — hrr— , Hi{ti) = - — — i j-,^^2(t2) = 1 + 2.1t2 and H^it^) = , ^ . In 

izt]^ 1 — zt ^ ~-t J-tS 

if we define A{ti) = iifi then ^ is a seventh-order iterative scheme (ALl) which is developed in |17| and 
if A{t\) = J- — then again 1 133 I I is a seventh-order iterative scheme (FD2-M1). 



IMPROVING THREE-POINT ITERATIVE METHODS FOR SOLVING NONLINEAR EQUATIONS 

Table 4. Numerical comparison between absolute errorsd^ri — a\^TNFE = 12) and COC 
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fnix),XO 


FD1-M2 


FD2-M1 


ALl 


FD1-M2 


FD2-M1 


ALl 










(COC) 


(COC) 


(COC) 


/i, 0.25 


3.60e-182 


8.55e-177 


6.80e-167 


6.9998 


6.9999 


6.9997 


/2, 1.1 


2.99e-263 


6.71e-172 


1.63e-145 


7.0000 


7.0000 


6.9999 


/3, 2.1 


8.60e-223 


3. 95c- 177 


3. lie- 178 


7.0000 


7.0000 


7.0000 


U, -0.5 


2.78e-297 


4.95e-247 


1.45e-233 


7.0000 


7.0000 


7.0000 


/s, 0.25 


1.31e-307 


1.45e-306 


2.81e-273 


7.0000 


7.0000 


7.0000 


/6, 1.2 


5.02e-304 


2.70e-246 


1.27e-222 


7.0000 


7.0000 


7.0000 


/7, 1.0 


1.51e-156 


1.83e-145 


3.98e-136 


6.9997 


6.9998 


6.9996 


/s, 1.6 


4.95e-240 


3.58e-102 


2.35e-67 


6.9995 


7.0001 


6.9978 


/9, 4.4 


5.10e-662 


3.74e-567 


1.07e-561 


7.0000 


7.0000 


7.0000 


/lO, 1-5 


1.08C-77 


1.33e-82 


1.73e-64 


6.8846 


7.0083 


6.9701 


/ii, 0.25 


8.67e-162 


4.51C-141 


1.31e-103 


7.0000 


7.0003 


7.0019 


/i2, 2.0 


1.30e-242 


4.29e-196 


9.45e-193 


7.0000 


7.0000 


6.9999 



Table 4 shows absolute error due to method FD2-M1 has small magnitude in comparison with other describe 
methods for the set of given twelve test functions and initial guesses. Only in the case of /g and fig, AL has better 
absolute error magnitude. 

3.3. Sixth-order convergence performance evaluation of FD4, FD5, FD6. For comparison, we develop 
some particular cases of derivative-free methods. In all calculations the value of parameter k = 1/100 and ti — ^^-'^"^ 



t2 



• FD4 



, is 



(65) 



(66) 



(67) 



^ Vn 



l+ti-t2 f(y,^) 



1-2*2 

2 



1 



-2t3-2t4 f[z„ 



1-2*2 1 
1-3*2 fU^,x„ 



• FD5 



Wn 


— 


l^f(Xn), 


Vn 






— 






= yn - 


1 


Zn 


l-*2 f[Xn,yrz\ ' 


Xn+1 




(, *5\" fM 



FD6 





— 


y-n 






= yn - 




— 



i-*i-H2 f(y^) 

1 — 2*1 fl=CTi,uJ„ 



• SK2M1, SK2M2 

In 1181 , Authors developed two methods by selecting different values of parameters. If a = /3 
then SK2 is SK2M1 and if o = /3 = r; = 1 then SK2 is SK2M2. 



Table 5 shows overall performance of derivative-free sixth-order iterative methods FD4, FD5, FD6. Absolute 
error \x„ — a\, for developed sixth-order derivative-free methods, is comparatively better than referenced sixth-order 
derivative-free iterative methods. Computational order of convergence is given in Table 6. dgt stands for divergent 
and X is for no information. 



3.4. Seventh-order convergence performance evaluation of FD7. FD7 is seventh-order convergent iterative 
method which is deduced from FD5 and is defined as 



(68) 





— 


Vn 


— Xji 


Zn 


= Vn 


Xn+1 


— Zfi 



■ l^fiXn), 

1-2*1 f(y^) 

1 — 3*1 flyn,11'n 
1 



l-t3 flZn.,V,t]+flZn..^n]-flx„,y„]- 
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Table 5. Numerical comparison between absolute errors(|a;n — a\,TNFE = 12) 



fnix) 


FD4 


FD5 


FD6 


TSl 


TS2 


SK2M1 


SK2M2 


FSl 


FS2 


fi, 0.25 


4.e-172 


2.e-158 


8.e-172 


9.e-54 


2.6-66 


3.6-103 


4.6-83 


1.6-59 


2.6-104 


/2, 1-1 


l.e-190 


2.e-297 


6.e-221 


0.01 


3.e-14 


7.6-162 


7.6-241 


2.e-7 


g.e-333 


/3, 2.1 


2.e-107 


2.e-203 


4.e-168 


3.e-18 


4.6-24 


2.6-197 


6.6-145 


3.6-7 


5.e-204 


/4, 0.5 


3.e-220 


2.e-198 


8.e-213 


3.e-127 


2.6-202 


1.6-217 


3.6-175 


9.6-196 


4.6-171 


/5, 0.25 


3.e-247 


7.e-248 


6.e-215 


2.e-129 


8.6-172 


4.6-157 


1.6-197 


2.6-171 


3.6-215 


/6, 1.2 


l.e-200 


2.e-186 


5.e-247 


l.e-132 


5.6-153 


8.6-112 


2.6-161 


8.6-194 


2.6-161 


/7, 1.0 


3.C-140 


8.C-130 


8.e-146 


dgt 


dgt 


3.6-78 


1.6-65 


dgt 


6.C-87 


/s, 1.6 


3.e-168 


l.c-90 


8.e-165 


5.C-55 


2.C-71 


2. 


5.6-25 


2.C-68 


3.C-145 


/9, 4.4 


5.e-466 


3.6-392 


1.6-404 


2.6-395 


2.C-424 


l.c-383 


l.c-366 


2.C-429 


5.6-378 


/lo, 1.5 


l.e-70 


6.6-70 


8.6-65 


l.e-138 


3.e-168 


0.06 


1. 


9.6-156 


5.6-68 


/ii, 0.25 


7.e-118 


4.6-118 


4.e-169 


2.6-71 


1.6-105 


1. 


8.6-29 


1.6-118 


3.C-85 


/l2, 2.0 


l.e-189 


3.e-149 


l.e-144 


8.6-137 


2.6-194 


2.6-173 


2.6-144 


5.e-215 


4.6-191 



Table 6. Computational order of convergence (COC) 





FD4 


FD5 


FD6 


TSl 


TS2 


SK2M1 


SK2M2 


FSl 


FS2 


h 


6.01 


6.00 


6.0 


5.99 


6.00 


5.98 


6.00 


5.98 


6.00 


h 


6.00 


6.00 


6.0 


-1.15 


6.24 


6.00 


6.00 


3.33 


6.00 


h 


6.00 


6.00 


6.0 


5.32 


5.63 


6.00 


6.00 


3.81 


6.00 


u 


6.00 


6.00 


6.0 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


h 


6.00 


6.00 


6.0 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


h 


6.00 


6.00 


6.0 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


h 


6.00 


6.00 


6.0 


X 


X 


5.98 


6.00 


X 


6.00 


h 


6.01 


6.00 


6.0 


5.99 


6.00 


5.23 


5.78 


6.00 


6.00 


h 


6.00 


6.00 


6.0 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 


fio 


6.30 


6.01 


6.0 


6.00 


6.00 


1.51 


0.00299 


6.00 


5.99 


fll 


6.05 


6.00 


6.0 


6.00 


6.00 


2.57 


5.91 


6.00 


6.00 


fl2 


5.98 


6.00 


6.0 


6.00 


6.00 


6.00 


6.00 


6.00 


6.00 



Table 7. Numerical comparison between absolute errors(|a;n — a\,TNFE = 12) and COC 





FD7(k)(COC) 


FS3-l(COC) 


FS3-2(COC) 


FS4-l(COC) 


FS4-2(COC) 


fl, 


0.25 


2.45e-378(1.0)(10) 


3.236-99(7) 


4.836-318(9) 


2.966-85(7) 


6.606-311(10) 


f2, 


1.1 


1.22e-388(0.01)(7) 


3.336-31(6.42) 


8.146-24(7.02) 


2.376-7(3.30) 


2.04(-0.0211) 


/3, 


2.1 


3.29e-266{0.01)(7) 


1.376-40(6.80) 


1.32c-26(7.42) 


3.96c-8(4.18) 


divergent (X) 


fi, 0.5 


8.61e-262(0.01)(7) 


1.896-221(7) 


1.17c-144(7) 


2.06c-210(7) 


8.75c-75(6.97) 


/5, 


0.25 


1.54c-303(0.01)(7) 


2.476-220(7) 


6.43e-397(7) 


1.336-206(7) 


1.256-384(7) 


/e, 


1.2 


5.15e-353(0.01)(7) 


1.71e-251(7) 


6.176-74(7.01) 


1.346-257(7) 


4.816-44(7.05) 


/7, 


1.0 


3.65e-213(0.01)(7) 


0.320(4.96) 


3.03(X) 


dgt(X) 


3.03(X) 


fa, 


1.6 


1.10e-174(0.01){7) 


1.09e-234(7) 


7.466-22(7.09) 


2.686-233(7) 


290(1.43) 


h, 


4.4 


9.O9e-689(0.01){7) 


4.556-617(7) 


1.716-517(7) 


5.066-640(7) 


4.456-506(7) 


fio, 


1.5 


3.81e-262(-1.0)(7) 


6.266-174(7) 


3.876-5(2.04) 


1.966-214(7) 


7.16(X) 


fll, 


0.25 


7.21e-212(0.01)(7) 


1.336-125(7) 


5.576-77(6.99) 


6.286-147(7) 


1.586-52(6.97) 


/l2, 


2.0 


3.52e-271(0.01)(7) 


7.936-257(7) 


6.79e6(-7.79) 


3.406-222(7) 


2.64(-0.261) 



All parameters set to zero value namely 7 = 0, 5 = in FS3-1, p = 0, t = in FS3-2 and w = 0, </> = in FS4. 
Table 7 represents the comparison between absolute errors as well as computational order of convergence. Clearly 
FD7 is superior than other methods in comparison. 

4. Conclusions 

We have presented derivative-based and derivative-free iterative methods by introducing free weighting parame- 
ters go, gi, 92 and weight functions. Weighting parameters, in the second step of derivative-free iterative methods, 
give flexibility for selecting different forms of approximation for flrst order derivative. A combination of derivative 
approximations can be constructed by fixing weighting parameters. The proper selection weight functions for both 
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derivative- free and derivative- based iterative methods help to improve the accuracy. The presented iterative meth- 
ods are defined in very general form which actually reproduce many existing iterative methods. By introducing 
different form of weighting functions and selection of parameters produce new families of iterative methods. We 
explore some of them to show the effectiveness of newly constructed iterative methods. 

Acknowledgement. This research was supported by Spanish MICINN grants AYA2010-15685 and MEC grants 
AYA2010-15685, AYA2008-04211-C02-C01. 
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